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Abstract 

We describe our recent results on the resonant perturbation theory of deco- 
herence and relaxation for quantum system with many qubits. The approach 
represents a rigorous analysis of the phenomenon of decoherence and relaxation 
for general TV-level systems coupled to reservoirs of the bosonic fields. We derive 
a representation of the reduced dynamics valid for all times t > and for small 
but fixed interaction strength. Our approach does not involve master equation 
approximations and applies to a wide variety of systems which are not explicitly 
solvable. 



1 Introduction 

Quantum computers (QCs) with large number of quantum bits (qubits) promise to 
solve important problems such as factorization of larger integer numbers, searching 
large unsorted databases, and simulations of physical processes exponentially faster 
than digital computers. Recently, many efforts were made for designing scalable (in 
the number of qubits) QC architectures based on solid-state implementations. One of 
the most promising designs of a solid-state QC is based on superconducting devices with 
Josephson junctions and solid-state quantum interference devices (SQUIDs) serving as 
qubits (effective spins), which operate in the quantum regime: frw » k^T, where T 
is the temperature and to is the qubit transition frequency. This condition is widely 
used in a superconducting quantum computation and quantum measurement, when 
T ~ 10- 20mK and Hu ~ 100- 150mK (in temperature units) [4,8,14-16,25-27] (see 
also references therein). The main advantages of a QC with superconducting qubits are: 
(i) The two basic states of a qubit are represented by the states of a superconducting 
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charge or current in the macroscopic (several [im size) device. The relatively large scale 
of this device facilitates the manufacturing, and potential controlling and measuring of 
the states of qubits. (ii) The states of charge- and current-based qubits can be measured 
using rapidly developing technologies, such as a single electron transistor, effective 
resonant oscillators and micro-cavities with RF amplifiers, and quantum tunneling 
effects, (iii) The quantum logic operation can be implemented exclusively by switching 
locally on and off voltages on controlling micro-contacts and magnetic fluxes, (iv) 
The devices based on superconducting qubits can potentially achieve large quantum 
dephasing and relaxation times of milliseconds and more at low temperatures, allowing 
quantum coherent computation for long enough times. In spite of significant progress, 
current devices with superconducting qubits only have one or two qubits operating 
with low fidelity even for simplest operations. 

One of the main problems which must be resolved in order to build a scalable QC 
is to develop novel approaches for suppression of unwanted effects such as decoherence 
and noise. This also requires to develop the rigorous mathematical tools for analyzing 
the dynamics of decoherence, entanglement and thermalization in order to control the 
quantum protocols with needed fidelity. These theoretical approaches must work for 
long enough times and be applicable to both solvable and not explicitly solvable (non- 
integrable) systems. 

Here we present a review of our results [18-20] on the rigorous analysis of the phe- 
nomenon of decoherence and relaxation for general iV-level systems coupled to reser- 
voirs. The latter are described by the bosonic fields. We suggest a new approach which 
applies to a wide variety of systems which are not explicitly solvable. We analyze in 
detail the dynamics of an iV-qubit quantum register collectively coupled to a thermal 
environment. Each spin experiences the same environment interaction, consisting of 
an energy conserving and an energy exchange part. We find the decay rates of the 
reduced density matrix elements in the energy basis. We show that the fastest decay 
rates of off-diagonal matrix elements induced by the energy conserving interaction is of 
order N 2 , while the one induced by the energy exchange interaction is of the order N 
only. Moreover, the diagonal matrix elements approach their limiting values at a rate 
independent of N. Our method is based on a dynamical quantum resonance theory 
valid for small, fixed values of the couplings, and uniformly in time for t > 0. We do 
not make Markov-, Born- or weak coupling (van Hove) approximations. 

2 Presentation of results 

We consider an TV-level quantum system S interacting with a heat bath R. The former 
is described by a Hilbert space f) s = and a Hamiltonian 



The environment R is modelled by a bosonic thermal reservoir with Hamiltonian 



H s = diag(£i,... 



E N ). 



(2.1) 




(2.2) 
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acting on the reservoir Hilbert space f)R, and where a*(k) and a(k) are the usual 
bosonic creation and annihilation operators satisfying the canonical commutation rela- 
tions [a(k), a* (I)] = S(k — I). It is understood that we consider R in the thermodynamic 
limit of infinite volume (1R 3 ) and fixed temperature T = l//3>0(ina phase without 
condensate). Given a form factor f(k), a square integrable function of k £ M 3 (momen- 
tum representation), the smoothed-out creation and annihilation operators are defined 
as a*(f) = f R3 f(k)a*(k)d 3 k and a(f) = J R3 f(k)a(k)d 3 k respectively, and the hermi- 
tian field operator is 

0(/) = -L [«*(/) + a(/)]. (2.3) 

The total Hamiltonian, acting on f)g ® F)r, has the form 

H = H s ® 1 R + Is ® H R + Xv, (2.4) 

where A is a coupling constant and v is an interaction operator linear in field operators. 
For simplicity of exposition we consider here initial states where S and R are not 
entangled, and where R is in thermal equilibrium. 1 The initial density matrix is thus 

Po = Po ® PR,/?, 

where p is any state of S and p^p is the equilibrium state of R at temperature 1/(3. 

Let A be an arbitrary observable of the system (an operator on the system Hilbert 
space f)s) and set 

(A) t := Tx s (p t A) = TT S+R (p t (A 1r)), (2.5) 
where pt is the density matrix of S + R at time t and 

Pt = TrRPt 

is the reduced density matrix of S. In our approach, the dynamics of the reduced 
density matrix ~p t is expressed in terms of the resonance structure of the system. Under 
the non- interacting dynamics (A = 0), the evolution of the reduced density matrix 
elements of S, expressed in the energy basis {tpk\k=i °f Hs, is given by 

[Ptlkl = ( V k,e- itHs p ^ tHs <Pi) = e iteifc [p ] kl , (2.6) 

where e\^ = E\ — E^. As the interaction with the reservoir is turned on, the dynamics 
(2.6) undergoes two qualitative changes. 

1. The "Bohr frequencies" 

e G {E — E' :E,E'e spec(F s )} (2.7) 

in the exponent of (2.6) become complex, e t— > e e . It can be shown generally 
that the resonance energies e e have non-negative imaginary parts, Ime e > 0. If 
Ime e > then the corresponding dynamical process is irreversible. 

x Our method applies also to initially entangled states, and arbitrary initial states of R normal w.r.t. 
the equilibrium state, see [19]. 
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2. The matrix elements do not evolve independently any more. Indeed, the effective 
energy of S is changed due to the interaction with the reservoirs, leading to a 
dynamics that does not leave eigenstates of He, invariant. (However, to lowest 
order in the interaction, the eigenspaces of H$ are left invariant and therefore 
matrix elements with (m, n) belonging to a fixed energy difference E m — E n will 
evolve in a coupled manner.) 

Our goal is to derive these two effects from the microscopic (hamiltonian) model and 
to quantify them. Our analysis yields the thermalization and decoherence times of 
quantum registers. 

2.1 Evolution of reduced dynamics of an iV-level system 

Let A € jB(I)s) be an observable of the system S. We show in [18, 19] that the ergodic 
averages 

((A)) ^ lim I [ T (A) t dt (2.8) 

exist, i.e., that (A) t converges in the ergodic sense as t — > oo. Furthermore, we show 
that for any t > and for any < uj' < 2ir/(3, 

(A) t ~ ((A)) „ = Y, j t£ Re(A) + O (aV^*) , (2.9) 

where the complex numbers e are the eigenvalues of a certain explicitly given operator 
K{ui'), lying in the strip {z £ C | < Imz < lo' /2}. They have the expansions 

e = e« = e + A 2 <^+0(A 4 ), (2.10) 

where e G spec(Hg (g) Is — Is <8> i?s) = s pec(i/s) — spec(ffs) and the 5^ are the 
eigenvalues of a matrix A e , called a level-shift operator, acting on the eigenspace of 
H$ Is — Is <8> -ffs corresponding to the eigenvalue e (which is a subspace of f)s (8> hs). 
The i? e (A) in (2.9) are linear functionals of A and are given in terms of the initial state, 
po, and certain operators depending on the Hamiltonian H. They have the expansion 

R £ (A)= X m ,nA m , n + 0(\ 2 ), (2.11) 

(m,n)6ie 

where I e is the collection of all pairs of indices such that e = E m — E n , the E^ being 
the eigenvalues of Hs- Here, A rritn is the (m,n)-matrix element of the observable A in 
the energy basis of Hs, and the x mi „ are coefficients depending on the initial state of 
the system (and on e, but not on A nor on A). 

Discussion. - In the absence of interaction (A = 0) we have e = e £ R, see 
(2.10). Depending on the interaction each resonance energy e may migrate into the 
upper complex plane, or it may stay on the real axis, as A / 0. 

- The averages (A) t approach their ergodic means ((A)) ^ if and only if Ime > for 
all e / 0. In this case the convergence takes place on the time scale [Ime] -1 . Otherwise 
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(A) t oscillates. A sufficient condition for decay is that Im<5e > (and A small, see 
(2.10)). 

- The error term in (2.9) is small in A, uniformly in t > 0, and it decays in time 
quicker than any of the main terms in the sum on the r.h.s.: indeed, Ime = 0(A 2 ) while 
u/ — 0(A) > lj'/2 independent of small values of A. However, this means that we are 
in the regime A 2 «. to' < 2ir/P (see before (2.9)), which implies that A 2 must be much 
smaller than the temperature T = 1/(3. Using a more refined analysis one can get rid 
of this condition, see also remarks on p. 376 of [19]. 

- Relation (2.11) shows that to lowest order in the perturbation the group of (energy 
basis) matrix elements of any observable A corresponding to a fixed energy difference 
E m — E n evolve jointly, while those of different such groups evolve independently. 

It is well known that there are two kinds of processes which drive decay (or irre- 
versibility) of S: energy-exchange processes characterized by [v, Hs] / and energy 
preserving ones where [v, Hs] = 0. The former are induced by interactions having 
nonvanishing probabilities for processes of absorption and emission of field quanta with 
energies corresponding to the Bohr frequencies of S and thus typically inducing thermal- 
ization of S. Energy preserving interactions suppress such processes, allowing only for a 
phase change of the system during the evolution ("phase damping", [3,7,10,13,22-24]). 

To our knowledge, energy-exchange systems have so far been treated using Born 
and Markov master equation approximations (Lindblad form of dynamics) or they have 
been studied numerically, while for energy conserving systems one often can find an 
exact solution. The present representation (2.9) gives a detailed picture of the dynamics 
of averages of observables for interactions with and without energy exchange. The 
resonance energies e and the functionals R £ can be calculated for concrete models, as 
illustrated in the next section. We mention that the resonance dynamics representation 
can be used to study the dynamics of entanglement of qubits coupled to local and 
collective reservoirs. Work on this is in progress. 

Contrast with weak coupling approximation. Our representation (2.9) of 
the true dynamics of S relies only on the smallness of the coupling parameter A, and 
no approximation is made. In the absence of an exact solution, it is common to make 
a weak coupling Lindblad master equation approximation of the dynamics, in which 
the reduced density matrix evolves according to p t = e t-C po> where C is the Lindblad 
generator, [2,5,6]. This approximation can lead to results that differ qualitatively 
from the true behaviour. For instance, the Lindblad master equation predicts that 
the system S approaches its Gibbs state at the temperature of the reservoir in the 
limit of large times. However, it is clear that in reality, the coupled system S + R will 
approach equilibrium, and hence the asymptotic state of S alone, being the reduction 
of the coupled equilibrium state, is the Gibbs state of S only to first approximation 
in the coupling (see also illustration below, and references [18,19]). In particular, the 
system's asymptotic density matrix is not diagonal in the original energy basis, but it 
has off-diagonal matrix elements of 0(A 2 ). Features of this kind cannot be captured 
by the Lindblad approximation, but are captured in our approach. 

It has been shown (see e.g. [5,6,9, 12]) that the weak coupling limit dynamics 
generated by the Lindblad operator is obtained in the regime A — > 0, t — ► oo, with 
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t/X 2 fixed. One of the strengths of our approach is that we do not impose any relation 
between A and t, and our results are valid for all times t > 0, provided A is small. 
It has been observed [9, 12] that for certain systems of the type S + R, the second 
order contribution of the exponents in (2.10) correspond to eigenvalues of the 
Lindblad generator. Our resonance method gives the true exponents, i.e., we do not 
lose the contributions of any order in the interaction. If the energy spectrum of H$ 

(s) 

is degenerate, it happens that the second order contributions to Imee vanish. This 
corresponds to a Lindblad generator having several real eigenvalues. In this situation 
the correct dynamics (approach to a final state) can be captured only by taking into 
account higher order contributions to the exponents ei s \ see [17]. To our knowledge, 
so far this can only be done with the method presented in this paper, and is beyond 
the reach of the weak coupling method. 

Illustration: single qubit. Consider S to be a single spin 1/2 with energy gap 
A = £2 — E\ > 0. S is coupled to the heat bath R via the operator 



a c 
c b 

where <p(g) is the Bose field operator (2.3), smeared out with a coupling function (form 
factor) g(k), k £ M 3 , and the 2x2 coupling matrix (representing the coupling operator 
in the energy eigenbasis) is hermitian. The operator (2.12) - or a sum of such terms, for 
which our technique works equally well - is the most general coupling which is linear 
in field operators. We refer to [19] for a discussion of the link between (2.12) and the 
spin-boson model. We take S initially in a coherent superposition in the energy basis, 



Po 



1 1 
1 1 



(2.13) 



In [19] we derive from representation (2.9) the following expressions for the dynamics 
of matrix elements, for all t > 0: 



[Pt\m,m 
[Pt\l,2 



-PEr, 



^ +^^tanh^^e ite °( A ) +J R m , m (A,t), m = 1,2, (2.14) 



l p ite-A(A) 



+ Rl,2{\t), 



(2.15) 



where the resonance energies e are given by 

eo(A) 
£a(A) 
£-a(A) 



®<Ks), (2-12) 



= iAV|c| 2 £(A) + 0(A 4 ) 

= A + \ 2 R+i\\ 2 [| c | 2 £(A) + (6-a) 2 £(0)] + 0(A 4 ) (2.16) 

= -^aTa) 

= lim- / d 3 kcoth(^] \g(k)\ 2 — ^ 5 , (2.17) 

40 ttJ r3 V 2 J m n {\k\ -7?) 2 + e 2 ' V ; 

R=\(b 2 -a 2 ){g,u- 1 g) 



with 
and 
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+±|c| 2 P.V. / u 2 | 5 (H,a)| 2 coth(M) _J_ . 

The remainder terms in (2.14), (2.15) satisfy |i? mi „(A, t)| < CA 2 , uniformly in t > 0, 
and they can be decomposed into a sum of a constant and a decaying part, 

e -HE m 

Rm,n{\ t) = {{p n>m )) - 5 mjn — h Bj {X, t), 

Z S,f3 

where \R' m n (\, t)\ = 0(A 2 e~ 7 *), with 7 = minjlmeo, Iul£±a}- These relations show 
that 

- To second order in A, convergence of the populations to the equilibrium values 
(Gibbs law), and decoherence occur exponentially fast, with rates tt = [Imeo(A)] -1 and 
td = [ImeA(A)] -1 , respectively. (If either of these imaginary parts vanishes then the 
corresponding process does not take place, of course.) In particular, coherence of the 
initial state stays preserved on time scales of the order A _2 [|c| 2 £(A) + (b — a) 2 £(0)] _1 , 
c.f. (2.16). 

- The final density matrix of the spin is not the Gibbs state of the qubit, and it is 
not diagonal in the energy basis. The deviation of the final state from the Gibbs state is 
given by lim^oo i? mi „(A,i) = 0(A 2 ). This is clear heuristically too, since typically the 
entire system S + R approaches its joint equilibrium in which S and R are entangled. 
The reduction of this state to S is the Gibbs state of S modulo 0(A 2 ) terms representing 
a shift in the effective energy of S due to the interaction with the bath. In this sense, 
coherence in the energy basis of S is created by thermalization. We have quantified 
this in [19], Theorem 3.3. 

- In a markovian master equation approach the above phenomenon (i.e., variations 
of 0(A 2 ) in the time-asymptotic limit) cannot be detected. Indeed in that approach 
one would conclude that S approaches its Gibbs state as t — > 00. 



2.2 Collective decoherence of a qubit register 

In the sequel we analyze in more detail the evolution of a qubit register of size N. The 
Hamiltonian is 

N N 

h s = j2JijSis*+J2 B i s h ( 2 - 18 ) 

i,j=l j=l 

where the Jij are pair interaction constants and Bj is the value of a magnetic field at 
the location of spin j. The Pauli spin operator is 



S z = 



1 

-1 



(2.19) 



and Sj is the matrix S z acting on the j-th spin only. 

We consider a collective coupling between the register S and the reservoir R: the 
distance between the iV qubits is much smaller than the correlation length of the 
reservoir and as a consequence, each qubit feels the same interaction with the reservoir. 
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The corresponding interaction operator is (compare with (2.4)) 

N N 

Aiui + \ 2 v 2 = Ai Y S- ® 4>{gi) + A 2 E ^(52). 



(2.20) 



Here 51 and 52 are form factors and the coupling constants Ai and A2 measure the 
strengths of the energy conserving (position-position) coupling, and the energy ex- 
change (spin flip) coupling, respectively. Spin-flips are implemented by the SJ in (2.20), 
representing the Pauli matrix 

" 1 



s x = 



1 



(2.21) 



acting on the j-ih spin. The total Hamiltonian takes the form (2.4) with Xv replaced 
by (2.20). It is convenient to represent p t as a matrix in the energy basis, consisting of 
eigenvectors ip^ of Hs- These are vectors in f)s = C 2 <8> • • • (g> C 2 = C 2 indexed by spin 
configurations 



a = {*!,...,**} €{+l,-l} N , 



fa = fu-i 



fcr N , 



where 



so that 





' 1 " 




' " 


(f + = 





. <p- = 


1 



N 



N 



H S ifa = E(a)(p s 



with 



We denote the reduced density matrix elements as 
The Bohr frequencies (2.7) are now 

N 



N 



e(a, t) = E{a) - E(t) = Y J ij( a i a j ~ T i T j) + Y B ^ a i ~ r i) 

i,j=l j=l 



(2.22) 
(2.23) 

(2.24) 

(2.25) 

(2.26) 



and they become complex resonance energies e e = £ e (Ai, A 2 ) € C under perturbation. 

Assumption of non- overlapping resonances. The Bohr frequencies (2.26) represent 
"unperturbed" energy levels and we follow their motion under perturbation (Ai,A 2 ). 
In this work, we consider the regime of non-overlapping resonances, meaning that the 
interaction is small relative to the spacing of the Bohr frequencies. 

We show in [19], Theorem 2.1, that for all t > 0, 



[Pt]a,r - «[PocW» = E eifee [ E W lr;a',r' \Po]*>*> + 0{\\ + A 

{e: e e ^0} o' ,r' 

+0((A 2 + A 2 )e-^'+°( A )] t ). 



(2.27) 
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This result is obtained by specializing (2.9) to the specific system at hand and con- 
sidering observables A = \(p T )((p a \. In (2.27), we have in accordance with (2.8) 

({\Poo\ ct,t)) = liniT-+oo t fo\Pt]a,T di. The coefficients w are overlaps of resonance 
eigenstates which vanish unless e = — e(a,r) = —e(a',T') (see point 2. after (2.7)). 
They represent the dominant contribution to the functionals R £ in (2.9), see also (2.11). 
The e e have the expansion 

e e = e^ = e + S^ +0(Xf + Xi), (2.28) 
where the label s = 1, . . . , v(e) indexes the splitting of the eigenavlue e into u(e) distinct 

is) 

resonance energies. The lowest order corrections <5e satisfy 

5<f) = 0{\\ + Al). (2.29) 

They are the (complex) eigenvalues of an operator A e , called the level shift operator 
associated to e. This operator acts on the eigenspace of L$ associated to the eigenvalue 
e (a subspace of the qubit register Hilbert space; see [19, 20] for the formal definition of 
A e ). It governs the lowest order shift of eigenvalues under perturbation. One can see 

(s) 

by direct calculation that Im5e > 0. 

Discussion. - To lowest order in the perturbation, the group of reduced density 
matrix elements [p t ]a,T associated to a fixed e = e(a, t) evolve in a coupled way, while 
groups of matrix elements associated to different e evolve independently. 

- The density matrix elements of a given group mix and evolve in time according 

■ , (s) 

to the weight functions w and the exponentials e e . In the absence of interaction 
(Ai = A2 = 0) all the e e s ^ = e are real. As the interaction is switched on, the ei^ 
typically migrate into the upper complex plane, but they may stay on the real line 
(due to some symmetry or due to an 'inefficient coupling'). 

- The matrix elements [p t ] _ ?) x of a group e approach their ergodic means if and only 

(s) 

if all the nonzero Ee have strictly positive imaginary part. In this case the convergence 
takes place on a time scale of the order l/7 e , where 

7e = min |lm4 s) : 8 = 1,..., v(e) s.t. e e s) + o} (2.30) 

is the decay rate of the group associated to e. If an e e s ^ stays real then the matrix 
elements of the corresponding group oscillate in time. A sufficient condition for decay 

(s) 

of the group associated to e is 7 e > 0, i.e. Im5 e > for all s, and Ai, A2 small. 

Decoherence rates. We illustrate our results on decoherence rates for a qubit 
register with Jij = (the general case is treated in [20] ) . We consider generic magnetic 
fields defined as follows. For rij G {0, ±1, ±2}, j = 1, . . . , N, we have 

TV 

B i n 3 = n i = V J- ( 2 - 31 ) 

Condition (2.31) is satisfied generically in the sense that only for very special choices 
of Bj does it not hold (one such special choice is Bj = constant). For instance, if the 
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Bj are chosen to be independent, and uniformly random from an interval [B m i n , B max \, 
then (2.31) is satisfied with probability one. We show in [20], Theorem 2.3, that the 
decoherence rates (2.30) are given by 

7e = ( i Me) + XlV2[e) + " 12(e) ' 6 + ? ) + 0(X\ + (2.32) 
[ A 2 2/ , e = J 

Here, y\ is a contributions coming from the energy conserving interaction, yo and yi 
are due to the spin flip interaction. The term y\2 is due to both interactions and is 
of 0(\\ + A|). We give explicit expressions for yo, yi, j/2 and yu in [20], Section 2. 
For the present purpose, we limit ourselves to discussing the properties of the latter 
quantities. 

- Properties ofyi(e): yi(e) vanishes if either e is such that eo := Y^j=i( a j~ T j) = 0> 
or the infra-red behaviour of the coupling function g\ is too regular (in three 
dimensions g\ oc \k\ p with p > —1/2). Otherwise yi(e) > 0. Moreover, y\(e) is 
proportional to the temperature T. 

- Properties of yi(e): 2/2 (e) > if g2(2Bj,T,) 7^ for all Bj (form factor g2(k) = 
g2(\k\,'E) in spherical coordinates). For low temperatures T, y2(e) oc T, for high 
temperatures ?/2 (e) approaches a constant. 

- Properties of y±2(e): If either of Ai, A2 or eo vanish, or if gi is infra-red regular 
as mentioned above, then y 12(e) = 0. Otherwise y 12(e) > 0, in which case y 12(e) 
approaches constant values for both T — > 0, 00. 

- FuW decoherence: If 7 e > for all e / then all off-diagonal matrix elements 
approach their limiting values exponentially fast. In this case we say that full 
decoherence occurs. It follows from the above points that we have full decoherence 
if A2 7^ and g2(2Bj, S) 7^ for all j, and provided Ai, A2 are small enough (so 
that the remainder term in (2.32) is small). Note that if A2 = then matrix 
elements associated to energy differences e such that eo = will not decay on the 
time scale given by the second order in the perturbation (Xf). 

We point out that generically, S + R will reach a joint equilibrium as t — > 00, 
which means that the final reduced density matrix of S is its Gibbs state modulo 
a peturbation of the order of the interaction between S and R, see [18, 19]. Hence 
generically, the density matrix of S does not become diagonal in the energy basis 
as t — ► 00. 

- Properties of yo-' Vo depends on the energy exchange interaction only. This 
reflects the fact that for a purely energy conserving interaction, the populations 
are conserved [18,19,23]. If g2(2Bj, S) / for all j, then yo > (this is sometimes 
called the "Fermi Golden Rule Condition"). For small temperatures T, yo oc T, 
while yo approaches a finite limit as T — > 00. 

In terms of complexity analysis, it is important to discuss the dependence of 7 e on 
the register size N. 
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- We show in [20] that y$ is independent of TV. This means that the thermalization 
time, or relaxation time of the diagonal matrix elements (corresponding to e = 0), 
is 0(1) in N. 

- To determine the order of magnitude of the decay rates of the off-diagonal density 
matrix elements (corresponding to e ^ 0) relative to the register size N, we 
assume the magnetic field to have a certain distribution denoted by ( ) . We show 
in [20] that 



(yi) = yi oc e\, (y 2 ) = C B ^{a - r), and (y 12 ) = c B (\i, A 2 )iV (e), (2.33) 



where Cb and cb = cb(Ai,A2) are positive constants (independent of N), with 
cb(Ai, A2) = 0{\\ + X 2 ). Here, iVo(e) is the number of indices j such that Uj = Tj 
for each {jj,t) s.t. e{a,r) = e, and 



is the Hamming distance between the spin configurations a and t (which depends 
on e only). 

- Consider e / 0. It follows from (2.32)-(2.34) that for purely energy conserving 
interactions (A2 = 0), 7 e oc Afeg = X 2 [Ylf=i(°~j ~ T j)] 2 i which can be as large as 
0(X 2 N 2 ). On the other hand, for purely energy exchanging interactions (Ai = 0), 
we have 7 e oc A|£>(_ct — r), which cannot exceed 0(X?,N). If both interactions are 
acting, then we have the additional term (yn), which is of order 0((X 2 + A|)iV). 
This shows the following: 

The fastest decay rate of reduced off- diagonal density matrix elements due to the 
energy conserving interaction alone is of order AfiV 2 , while the fastest decay rate 
due to the energy exchange interaction alone is of the order X^N. Moreover, the 
decay of the diagonal matrix elements is of oder \\, i.e., independent of N. 

Remarks. 1. For A2 = the model can be solved explicitly [23], and one shows that 
the fastest decaying matrix elements have decay rate proportional to X 2 N 2 . Further- 
more, the model with a non-collective, energy- conserving interaction, where each qubit 
is coupled to an independent reservoir, can also be solved explicitly [23]. The fastest 
decay rate in this case is shown to be proportional to Af N. 

2. As mentioned at the beginning of this section, we take the coupling constants 
Ai, A2 so small that the resonances do not overlap. Consequently X 2 N 2 and X 2 N are 
bounded above by a constant proportional to the gradient of the magnetic field in the 
present situation, see also [20]. Thus the decay rates 7 e do not increase indefinitely 
with increasing N in the regime considered here. Rather, the 7 e are attenuated by 
small coupling constants for large N. They are of the order 7 e ~ A. We have shown 
that modulo an overall, common (iV-dependent) prefactor, the decay rates originating 
from the energy conserving and exchanging interactions differ by a factor N. 



N 




(2.34) 
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3. Collective decoherence has been studied extensively in the literature. Among the 
many theoretical, numerical and experimental works we mention here only [1,3,10,11, 
23], which are closest to the present work. We are not aware of any prior work giving 
explicit decoherence rates of a register for not explicitly solvable models, and without 
making master equation technique approximations. 

3 Resonance representation of reduced dynamics 

The goal of this section is to give a precise statement of the core representation (2.9), 
and to outline the main ideas behind the proof of it. 

The iV-level system is coupled to the reservoir (see also (2.1), (2.2)) through the 
operator 

R 

v = ^ X r G r <8> 4>(g r ), (3.1) 

r=l 

where each G r is a hermitian N x N matrix, the g r {k) are form factors and the A r eR 
are coupling constants. Fix any phase x £ ^ an d define 

i \ I ^ I 1 1/2 f 9r(u, c) if u > 0, , n n . 

■= Vr^T I-I 1/2 { -eV^(-u,a) if u < o, ( 3 - 2 ) 

where u G M and a G S 2 . The phase % is a parameter which can be chosen appropriately 
as to satisfy the following condition. 

(A) The map oj i— > g r g{u + u, a) has an analytic extension to a complex neighbour- 
hood {\z\ < uj'} of the origin, as a map from C to L 2 (IR 3 , d 3 k). 

Examples of g satisfying (A) are given by g(r,a) = r p e _T * m <?i(<x), where p = —1/2 + 
n, n = 0, 1, . . ., m = 1, 2, and gi(cr) = e 1 ^g 1 (a). 

This condition ensures that the technically simplest version of the dynamical reso- 
nance theory, based on complex spectral translations, can be implemented. The tech- 
nical simplicity comes at a price: on the one hand, it limits the class of admissible 
functions g(k), which have to behave appropriately in the infra-red regime so that the 
parts of (3.2) fit nicely together at u = 0, to allow for an analytic continuation. On 
the other hand, the square root in (3.2) must be analytic as well, which implies the 
condition lo' < 2tt/(3. 

It is convenient to introduce the doubled Hilbert space Tts = f)s <8> f)s 5 whose nor- 
malized vectors accommodate any state on the system S (pure or mixed). The trace 
state, or infinite temperature state, is represented by the vector 

1 N 

ns = ~nu Yl W ® ¥>j ( 3 - 3 ) 

via 

B(t>s)3A^(n s ,{A®t)n s ). (3.4) 

Here the <pj are the orthonormal eigenvectors of H$. This is just the Gelfand-Naimark- 
Segal construction for the trace state. Similarly, let and fl^g be the Hilbert 
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space and the vector representing the equilibrium state of the reservoirs at inverse 
temperature (3. In the Araki- Woods representation of the field, we have Wr = T ® T , 
where T is the bosonic Fock space over the one-particle space L 2 (M 3 ,d 3 /c) and Or^ = 
Q, O being the Fock vacuum of T (see also [19,20] for more detail). Let ipo <8> £l-R t p 
be the vector in Tig ® Wr representing the density matrix at time t = 0. It is not 
difficult to construct the unique operator in B G Us <8> f)s satisfying 

Bn s = Vo. 

(See also [19] for concrete examples.) We define the reference vector 

O re f := f2 s ® ^R/3 

and set 

A = max I X r I . 

r=l,...,fl 

Theorem 3.1 (Dynamical resonance theory [18-20]) Assume condition (A) with 
a fixed uj' satisfying < to' < 2ir/f3. There is a constant cq s.t. if A < cq/P then the 
limit ((A)) (2.8), exists for all observables A G B(t)s)- Moreover, for all such A and 
for all t > we have 

<A-«^»oc = E E eite/ ((B*i>o)®n R ,p,Qi s \A®i s )n re i) 

e,s: ei s) ^0 8=1 

+ 0(X 2 e^' +0 ^ t ). (3.5) 

The ei s ^ are given by (2.10), 1 < v(e) < mult(e) counts the splitting of the eigenvalue 
e into distinct resonance energies e{ s ^ and the are (non- orthogonal) finite-rank 

projections. 

This result is the basis for a detailed analysis of the reduced dynamics of concrete 
systems, like the iV-qubit register introduced in Section 2.2. We obtain (2.27) (in 
particular, the overlap functions w) from (3.5) by analyzing the projections in 
more detail. Let us explain how to link the overlap (^(B*ipo) <£> ^r^, Q^\a ® is)^ re f^) 

to its initial value for a non-degenerate Bohr energy e, and where A = |v3 n )( ( ^ m |. (The 
latter observables used in (2.9) give the matrix elements of the reduced density matrix 
in the energy basis.) 

The is the spectral (Riesz) projection of an operator K\ associated with the 
eigenvalue e{ s \ see (3.11). 2 If a Bohr energy e, (2.7), is simple, then there is a single 

(s) 

resonance energy e e bifurcating out of e, as A / 0. In this case the projection Q e = Q e 
has rank one, Q e = |Xe)(Xe|> where \ e and \e are eigenvectors of K\ and its adjoint, 
with eigenvalue e e and its complex conjugate, respectively, and {Xe,Xe) = 1- From 

2 In reality, we consider a spectral deformation K\(u>), where u> is a complex parameter. This is a 
technical trick to perform our analysis. Physical quantities do not depend on w and therefore we do 
not display this parameter here. 
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perturbation theory we obtain Xe = Xe = <8> fi <8> ^r,/3 + 0(A), where -ffsVj = ^j^j 
and Ek — Ei = e. The overlap in the sum of (3.5) becomes 

(0B>o) ® ^R,/3, Qe(A ® l S Ref) 

= <(B*^o) ® fi Ri/3 , |p fc ® ® ^R/3>(^fc ® ¥?j (8) n Rlj 8|(4 ® IsRcf) + 0{\ 2 ) 
= < J B*^ ,|yfc®yj)<yfc®yj|(^®ls)ns) + O(A 2 ) (3.6) 

The choice A = |<^n)( ( ^ , m| in (2-5) gives ( J 4) i = [p t ] m , n , the reduced density matrix 
element. With this choice of A, the main term in (3.6) becomes (see also (3.3)) 

(B*ip , \ip k ® <pi)(<p k ® ® is)n s ) 

V -/V 

= 4n<W (-8*^0, (|¥>n)(¥>m| <8> ls)^S> 
= 4n^m (V>0, (|<Pn)(<£m| <8> ls)-BO S ) 

= <Jfcn^ m [po]mra- (3-7) 

In the second-last step, we commute B to the right through \(p n ){'-Pm\ <8> Us, since B 
belongs to the commutant of the algebra of observables of S. In the last step, we use 
B£l s = Vo- 

Combining (3.6) and (3.7) with Theorem 3.1 we obtain, in case e = E m — E n is a 
simple eigenvalue, 

[PtUn ~ «D>ooW» = E ^ [tkntlmtPolmn + 0(\ 2 )] + 0(AV^ +0 ^). 

This explains the form (2.27) for a simple Bohr energy e. The case of degenerate e 
(i.e., where several different pairs of indices k,l satisfy E k — E\ = e) is analyzed along 
the same lines, see [20] for details. 

3.1 Mechanism of dynamical resonance theory, outline of proof of 
Theorem 3.1 

Consider any observable A € -B(f)s)- We have 

(A) t = Tr s [p t A] 

= Trs+R [pt A ® 1 R ] 

= (^o,e itL M^®ls®lR]e- itL ^>. (3.8) 

In the last step, we pass to the representation Hilbert space of the system (the GNS 
Hilbert space), where the initial density matrix is represented by the vector i/jq (in 
particular, the Hilbert space of the small system becomes f)s <8> f)s)j see also before 
equations (3.3), (3.4). As mentioned above, in this review we consider initial states 
where S and R are not entangled. The initial state is represented by the product vector 
ipo = fi s ® ^R,/3, where £l s is the trace state of S, (3.4), (fi s , (^4 ® is)^s) = Ty-Tr (A), 
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and where ^R j( g is the equilibrium state of R at a fixed inverse temperature < [3 < oo. 
The dynamics is implemented by the group of automorphisms e ltLx • e~ ltLx . The self- 
adjoint generator L\ is called the Liouville operator. It is of the form L\ = Lq + XW, 
where Lq = Ls + represents the uncoupled Liouville operator, and XW is the 
interaction (3.1) represented in the GNS Hilbert space. We refer to [19, 20] for the 
specific form of W. 

We borrow a trick from the analysis of open systems far from equilibrium: there is 
a (non-self-adjoint) generator K\ s.t. 

e itLA Ae~ itLx = e itKx Ae~ itKx for all observables A,t>0, and 
Kxipo = 0. 

K\ can be constructed in a standard way, given L\ and the reference vector ip . K\ is 
of the form K\ = Lq + XI, where the interaction term undergoes a certain modification 
(W — ► I), c.f. [19]. As a consequence, formally, we may replace the propagators in 
(3.8) by those involving K. The resulting propagator which is directly applied to ipo 
will then just disappear due to the invariance of ipQ. One can carry out this procedure 
in a rigorous manner, obtaining the following resolvent representation [19] 

(A) t = ~ f (ik,(K x (u)-z)- 1 [A®l a ®l R ]fa)tP*dz, (3.9) 
2iri J R _i 

where K\(lu) = Lq(uj) + XI (uj), I is representing the interaction, and uj \— > K\{uj) is a 
spectral deformation (translation) of K\. The latter is constructed as follows. There 
is a deformation transformation U{ui) = e~ lujD , where D is the (explicit) self-adjoint 
generator of translations [19-21] transforming the operator K\ as 

K x (lo) = U^KxUiio)- 1 = L + loN + XI{lo). (3.10) 
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Figure 1: Spectrum of Kq(u>) 

Here, iV = N\ (8) 1+ l<g>Ni is the total number operator of a product of two bosonic 
Fock spaces T ® T (the Gelfand-Naimark-Segal Hilbert space of the reservoir), and 
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Figure 2: Spectrum of K\(ui). Resonances e£ are uncovered. 

where N\ is the usual number operator on T . N has spectrum N U {0}, where is 
a simple eigenvalue (with vacuum eigenvector Qr/j = f2 ® O). For real values of w, 
U(u) is a group of unitaries. The spectrum of K\(u>) depends on Imu and moves 
according to the value of Imu;, whence the name "spectral deformation". Even though 
U(u) becomes unbounded for complex u>, the r.h.s. of (3.10) is a well defined closed 
operator on a dense domain, analytic in uj at zero. Analyticity is used in the derivation 
of (3.9) and this is where the analyticity condition (A) after (3.2) comes into play. The 
operator I(u>) is infinitesimally small with respect to the number operator N. Hence 
we use perturbation theory in A to examine the spectrum of K\{oS). 

The point of the spectral deformation is that the (important part of the) spectrum 
of K\{oj) is much easier to analyze than that of K\, because the deformation uncovers 
the resonances of K\. We have (see Figure 1) 

spec (Ko(u)) = {Ei - /•.', },., i v |J {ra + R}, 

n>l 

because Kq{lj) = Lq + loN, Lq and N commute, and the eigenvectors of Lq = L§ + Lr 
are ifi (g> (fj (8> ^r,^- Here, we have H$(fj = Ejipj. The continuous spectrum is bounded 
away from the isolated eigenvalues by a gap of size Imw. For values of the coupling 
parameter A small compared to Imw, we can follow the displacements of the eigenvalues 
by using analytic perturbation theory. (Note that for Imu = 0, the eigenvalues are 
imbedded into the continuous spectrum, and analytic perturbation theory is not valid! 
The spectral deformation is indeed very useful!) 

Theorem 3.2 ( [19]) (See Fig. 2.) Fix Imw s.t. < Imw < oj' (where uj' is as 
in Condition (A)). There is a constant cq > s.t. if |A| < cq//3 then, for all u 
with Imw > 7w'/8, the spectrum of K\{uj) in the complex half-plane {Imz < w'/2} is 
independent of to and consists purely of the distinct eigenvalues 

: e G spec(L s ),s = 1, . . . ,i/(e)}, 

where 1 < v(e) < mult(e) counts the splitting of the eigenvalue e. Moreover, 

lim | e W(A)-e| = 
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for all s, and we have Im4 > 0. Also, the continuous spectrum of K\(u) lies in the 
region {Imz > 3u//4}. 

Next we separate the contributions to the path integral in (3.9) coming from the 
singularities at the resonance energies and from the continuous spectrum. We deform 
the path of integration z = R — i into the line z = R + iu//2, thereby picking up the 
residues of poles of the integrand at si ( & H e, s). Let be a small circle around ei\ 
not enclosing or touching any other spectrum of K\(uj). We introduce the generally 
non-orthogonal Riesz spectral projections 




(3.11) 



Figure 3: Contour deformation: J R _. dz = Y^ e , s I c i s) ^ z + fft+iu' /2 ^ z 
It follows from (3.9) that 

u(e) 

= E E ^ (^o, Qi s) [A®ts® IrJ^o) + 0(A 2 e-'*/ 2 ). (3.12) 

e s=l 

Note that the imaginary parts of all resonance energies si' are smaller than ui' /2, so 
that the remainder term in (3.12) is not only small in A, but it also decays faster than 
all of the terms in the sum. (See also Figure 3.) We point out also that instead of 
deforming the path integration contour as explained before (3.11), we could choose 
z = R + i[u/ — 0(A)], hence transforming the error term in (3.12) into the one given in 
(3.5). 
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Finally, we notice that all terms in (3.12) with ei / will vanish in the ergodic 
mean limit, so 



This concludes the outline of the proof of Theorem 3.1. 
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